Pair-density wave states through spin-orbit coupling in multilayer superconductors 
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Spin singlet superconductors with quasi-two dimensional multilayer structure are studied in high 
magnetic fields. Specifically we concentrate on bi- and tri-layer systems whose layers by symme- 
try are subject Rashba-type spin-orbit coupling. The combination of magnetic field and spin-orbit 
coupling leads to a first order phase transition between diflterent states of layer-dependent supercon- 
ducting order parameters upon rising the magnetic field. In this context we distinguish the low-field 
Bardeen-Cooper-Schrieffer state where all layers have order parameters of the same sign and the 
high-field pair-density wave state where the layer-dependent order parameters change the sign at the 
center layer. We also show that progressive paramagnetic limiting effects yield additional features 
in the H-T phase diagram. As possible realizations of such unusual superconducting phases we 
consider artificial super lattices of CeCoIns, as well as some multilayer high-Tc cuprates. 

PACS numbers; 74.20.-z, 74.25.Dw, 74.70.Tx 



I. INTRODUCTION 

Time reversal and inversion are two of the most im- 
portant symmetries for the formation of Cooper pairs 
in a superconductor. Interesting new features of super- 
conductivity arise when one or both of them are missing. 
One intriguing phase in the absence of time reversal sym- 
metry is the so-called Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) state - a spatially modulated superconducting 
condensate induced by the spin-splitting of the Fermi 
surface due to high magnetic fields;-'^ This phase has 
been discussed in many contexts, for condensed matter 
systems^i^ and ultra-cold atomic gases^ as well as in nu- 
clear physics The other interesting class is represented 
by superconductors without inversion center in the crys- 
tal lattice. These so-called non-centrosymmetric super- 
conductors are also characterized by spin-splitting of the 
Fermi surface, in this case due to anti-symmetric spin- 
orbit coupling, and form Cooper pairs of mixed parityJ 
Interestingly, both types of situations have been recently 
discussed in the context of Ce-based heavy Fermion su- 
perconductors. A phase observed in CeCoIns at low- 
temperatures and rather high magnetic fields is consid- 
ered a possible candidate for an FFLO state.— i^i^i^ Most 
striking properties of unusually high upper critical fields 
have been reported for CelrSis and CeRhSia which both 
have a non-centrosymmetric crystal lattice.^ In our study 
we would like to consider a combined effect of spin- 
splitting through magnetic fields and spin-orbit coupling 
in an artificially structured heavy Fermion superconduc- 
tor based on CeCoIus. 

Mizukami and his coworkers succeeded recently in fab- 
ricating superlattices of CeCoIns and YbCoIns, the for- 
mer a heavy Fermion and the latter an ordinary metal 
Surprisingly this system is superconducting, if the num- 
ber of CeCoIns layers within a unit cell exceeds two. 
In view of the fact that the upper critical field re- 
mains amazingly high despite the low transition tempera- 
ture^ the superconductivity in the superlattice is robust 



against paramagnetic depairing effect. It was proposed 
that the superstructure of this material would have lo- 
cally non-centrosymmetric features while having global 
inversion centers and may lead to unusual properties of 
the superconducting phase. Actually the multilayered 
structure would induce a spatially dependent Rashba- 
type of spin-orbit couplingi^ which would influence the 
spin structure of the electronic bands and suppress the 
paramagnetic depairing effect.^- Recent measurement of 
the angular variation of the upper critical field found evi- 
dence for such a behavior In this study, we investigate 
an unusual form of the superconducting order parameter 
induced by the modulated Rashba spin-orbit coupling in 
the magnetic field. As it will turn out to involve a mod- 
ulation on the length scale of lattice constant, the phase 
may be viewed as a pair-density wave (PDW) state. Here 
we will restrict to models of bi- and tri-layer systems as 
depicted schematically in Fig. [TJ 
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FIG. 1. (Color onfine) Schematic figures of (a) bi- and (b) tri- 
layer system. Thick bars describe the two dimensional con- 
ducting planes which are coupled with each other through the 
inter-layer coupling t±. The structures of modulated Rashba 
spin-orbit coupling and superconducting order parameters in 
BCS and PDW states are shown in the right hand side of each 
figure. See the text for details. 



II. FORMULATION 

For our discussion we analyze the following model of a 
multilayered system having the inter-layer hopping and 
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Rashba spin-orbit coupling, 

k^s.m fc,s,(m,m') 
fc,fc',m 

-(J-bH Sci^^Cksm, (1) 

k,s,m 

where m is the index of layers. Note that we restrict here 
to one set of superconducting layers and ignore the nor- 
mal metal part of the superlattice, assuming that the cou- 
pling between heavy Fermion and normal metal part is 
small. This is a reasonable approximation for the super- 
lattice of CeCoIns , since the large mismatch in the Fermi 
velocities between CeCoIns and YbCoIns at the interface 
suppresses the proximity effect .~ For the band structure 
we assume a simple nearest-neighbor hopping tight bind- 
ing form ^(fc) = — 2i(cosA;a; -|- cosky) — on a square 
lattice with the chemical potential fi/t — 2, and a small 
inter-layer hopping t±/t = 0.1. The symbol (m,m') de- 
notes the summation over nearest-neighbor layers. It is 
important to note that the following results do not quali- 
tatively depend on the details of the band structure. We 
also use an s-wave superconducting state by assuming 
V{k, k') = ~V (onsite attractive interaction) for sim- 
plicity. Although CeCoIns is believed to be a d-wave 
superconductor and would have a small odd-parity com- 
ponent admixed with the local non-centrosymmetricity, 
the features found below are independent of the pairing 
symmetry as long as spin-singlet pairing is predominant. 
In fact, we have confirmed that the model for d-wave su- 
perconductivity gives the qualitatively same results (see 
Appendix). Hence, for the sake of numerical accuracy, 
we restrict ourselves to the simpler s-wave state. Our 
choice V/t — 1.7 gives rise to the critical temperature 
kBTco/t = 0.0255 for both bi-layer and tri-layer systems 
with a = 0, ignoring phase fluctuations. 

For our purpose it is sufficient to choose a g- vector with 
Rashba spin-orbit structure, g{k) = (— sin fcy, sin fca;, 0) 
without going into microscopic details of the com- 
pound!^ The coupling constants am are layer dependent 
and antisymmetric with respect to reflection at the cen- 
ter of the multilayer structure. In this way the existence 
of a global inversion center is guaranteed. For exam- 
ple, (ai,a2) = (a, —a) for bi-layers, and (ai, 0:2,03) — 
(a, 0,-0;) for tri-layers (see Fig. [1]). Paramagnetic lim- 
iting and also eventually the PDW phase enter through 
the Zeeman coupling term, i.e. the last term in Eq. ([l}. 

In our discussion we restrict to magnetic fields 
along c-axis and analyze the model on the basis of 
the Bogoliubov-de Gennes equation. While we in- 
clude the layer dependence of order parameter = 
— J2k^i^-kimCkf7n), we ignore the spatial modulation 
in the plane, such as the vortex state. This simplification 



is not crucial on a qualitative level, if we assume short 
coherence lengths and a large Ginzburg-Landau param- 
eter K, the ratio of coherence length and London pene- 
tration depth. Although phase fluctuations may play an 
important role for the critical temperature, we focus on 
the low-temperature and high-magnetic field phase, and 
ignore this aspect. 



III. SUPERCONDUCTING PHASES 
A. Bi-layer system 

First, we address the results for the bi-layer system. 
Figure [5] shows the T-H phase diagram for several val- 
ues of a/t±. For a/t± — [Fig. UKa)], we obtain the 
conventional phase diagram of a superconductor subject 
to paramagnetic limiting^S, In this case the order param- 
eter is the same in both layers having the same phase, 
(Ai,A2) = (A, A). This state we call Bardeen-Cooper- 
Schrieffer (ECS) state is stabilized basically by the inter- 
layer "Josephson" coupling. The superconducting phase 
transition turns first order at low temperatures, because 
of the paramagnetic depairing effect. The in-plane FFLO 
phase with an inhomogeneous order parameter in the 
two-dimensional plane is supposed to appear around the 
first order transition line in Fig.[2Ka). However, in our ap- 
proach we ignore the possibility of in-plane modulations 
as mentioned above. 

We now turn to the role of the layer-dependent Rashba 
spin-orbit coupling. Interestingly, the order parameter 
configuration with a sign change, (Ai,A2) = (A, —A), 
becomes more stable at high magnetic fields ioi a/t± > 1. 
This state is called the pair-density wave (PDW) state 
in Figs, mb-d), since the order parameter modulates 
on length scales of the crystal lattice constant. The 
PDW state was proposed for the Q-phase in the bulk 
CeCoIus,— but has not been established in microscopic 
models so far. With increasing the spin-orbit coupling a 
the PDW state is stabilized in the large parameter range. 
Thus, we here flnd that the PDW state may be real- 
ized in multilayer systems having "weak inter-layer cou- 
pling" and "moderate spin-orbit coupling". We would 
like to stress that the PDW state is stabilized in the 
purely superconducting multilayers, while a similar in- 
homogeneous superconducting state can be induced by 
the proximity effect in the ferromagnet-superconductor 
junctionJ^ 

The stability of the PDW state can be understood on 
the basis of the spin susceptibility and the band structure 
investigated in Ref. [ill In the absence of spin-orbit cou- 
pling the bi-layer system forms two bands, the bonding 
band and anti-bonding band. Cooper pairs in the PDW 
state would form by inter-band pairing, in contrast to 
the BCS state which is an intra-band pairing state. Ob- 
viously, the PDW state would not be stable under these 
circumstances. On the other hand, the band structure is 
described by the two spin-split bands on each layer for 
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FIG. 2. T-H phase diagram of bi-layer systems for (a) a/t± = 0, (b) a/t±_ = 1, (c) a/t± — 2, and (d) a/t± = 3. The sohd and 
dashed hnes show the second- and first-order phase transition lines, respectively. The "BCS" and "PDW" are explained in the 
text. The temperature T and magnetic field hbH in Fig. [2]are normalized by the critical temperature Tco at a/t± = 0. 



large a/t±. In this case, both BCS and PDW states occur 
as intra-band Cooper pairing states. Although the zero- 
field transition temperature of the PDW state is lower 
than the onset of superconductivity (BCS state), a mag- 
netic field along the c-axis suppresses the BCS state more 
strongly through paramagnetic depairing. Indeed, the 
PDW state has a larger spin susceptibility than the BCS 
state.-ii Therefore, the PDW state is more robust against 
paramagnetic limiting than the BCS state, as shown in 
Figs.Hb-d). 



B. Tri- layer system 

New aspects can be found for the tri-layer system. 
The order parameter has now three components, such 
that we distinguish the BCS state as (Ai,A2,A3) — 
(Aout, Ain, Aout) and the PDW state as (Ai,A2, A3) = 
(Aout, 0, — Aout)- For a/t± < 2, the phase diagram is 
qualitatively the same as shown in Fig.[2]for bi-layer sys- 
tems. The range of the PDW state is slightly reduced 
because of the complete suppression of order parame- 
ter in the inner layer. Increasing spin-orbit coupling to 
a/t± = 3 yields an intriguing modification of the phase 
diagram (see Fig. [3]). An further first-order phase tran- 
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FIG. 3. T-H phase diagram of tri-layer systems for a/t± — 3. 
The dash-dotted line shows the crossover field determined by 
the criterion Ain = Aout- 

sition line appears within the BCS state, which is absent 
for bi-layer systems. 

In order to characterize the first order phase transition 
in the BCS state, we show the magnetic field dependence 
of the order parameter in Fig.|4l At H = 0, the order pa- 
rameter of the inner layer Ajn is slightly larger than that 
of the outer layers Aout- Upon increasing the magnetic 
field for r/Tco = 0-04 at usH/kBTco = 2-04 an abrupt 
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FIG. 4. (Color online) Magnetic field dependence of the order 
parameters in the inner layer Ain (solid lines) and outer layers 
Aout (dashed lines) for a/tx = 3. (a) T/Tco = 0.04 and (b) 
T/Tco = 0.28. 



drop of Ain occurs, while Aout is only weakly affected [see 
Fig. m^a)]. Indeed, with the layer-dependent spin-orbit 
coupling the paramagnetic limiting behavior is different 
for the inner and outer layers. The spin-susceptibility in 
the superconducting phase is more strongly diminished 
for the inner than the outer layer when a/t j_ > I , as the 
Rashba spin-orbit coupling vanishes by symmetry for the 
inner layeiJ^ The high-field BCS state is stabilized by 
the more field resistant order parameter component Aout 
of the outer layer. Note that the inter-layer hopping is 
responsible for the fact that the order parameter compo- 
nents are not behaving independently. As the tempera- 
ture is increased, the abrupt first order phase transition 
in the BCS state turns into a crossover at the critical 
end point T/Tco ~ 0.12. In this regime the order param- 
eter of the inner layer decreases continuously, as shown 
in Fig.mjb). 

The different phases also influence the quasiparticle 
spectrum which we will discuss here looking at the den- 
sity of states (DOS) layer- resolved for the tri-layer sys- 
tem, as depicted in Fig. [5] (solid line for inner and dashed 
line for outer layer). The DOS of the outer layer always 
shows a gap at the Fermi energy. In the absence of a 
magnetic field the DOS of both the inner and outer lay- 




FIG. 5. (Color online) DOS in the inner layer pin{oj) (solid 
lines) and outer layers /9out(i^) (dashed lines) at T/Tco = 0.04. 
(a) /iBH/ksTco = (low field BCS state), (b) usH/kBTco = 
2.25 (high field BCS state), and (c) fiBH/kBTco = 3 (PDW 
state), respectively. The other parameters are the same as 
Fig. [SI 



ers are essentially identical [see Fig. O^a)]. Differences 
appear for finite fields. In Fig. ^h) we show the DOS 
in the high-field BCS state, which is different for the in- 
ner and outer layers. The quasiparticle spectrum of the 
inner layer is gapless, although the order parameter Ajn 
does not disappear. Surprisingly, a small gap on the in- 
ner layer opens up again when we enter the PDW state 
as revealed by Fig. [5jc), although the order parameter 
Ain vanishes in this case. The origin of this somewhat 
counter-intuitive behavior lies in the inter-layer hopping 
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and the energy shifts due to the magnetic field leading to 
the pair breaking. The inner layer is obviously affected by 
paramagnetic limiting in the high-field BCS state yield- 
ing subgap quasiparticle states. On the other hand, in the 
PDW state the paramagnetic pair breaking effects almost 
disappear?^ In this case the intrinsic order parameter of 
the inner layer vanishes and a small gap appears in the 
quasiparticle excitations because the quasiparticles are 
extended over all layers due to the inter-layer hopping. 



IV. SUMARRY AND DISCUSSION 

As mentioned earlier, possible realizations of the situ- 
ation discussed here could be found in artificially grown 
superlattices of CeCoIng^S as well as in multilayer high- 
Tc cuprates.— 1^ These multilayer systems are quasi- two- 
dimensional superconductors and likely near the Pauli 
limit as they are rather robust against orbital depair- 
ing due to the very short coherence lengths . Conse- 
quently, they satisfy the condition for the phenomenol- 
ogy analyzed in our study. Indeed, the reduced param- 
agnetic limiting due to the Rashba spin-orbit coupling of 
the outer layers has been confirmed in recent experiments 
for superlattices CeCoIns,^- The experiment is approach- 
ing to the FFLO state for the field along a6-planeJ^ Our 
study shows that the PDW state may be stabilized in the 
magnetic field along c-axis. While this lies in the mea- 
surable range for the "low-Tc" system CeCoIus, fields 
revealing the high-field BCS or PDW state are likely be- 
yond 100 T in high-Tc cuprates. 

Finally, we would like to comment on several points of 
our results. First, we restricted to pairing interaction in 
the spin singlet channel. Spin-singlet pairing is essential 
for the observation of the phase diagram including differ- 
ent BCS states and the PDW state. If spin triplet pairing 
compatible with the Rashba-spin-orbit coupling would be 
dominant, paramagnetic limiting would not play a role at 
all, since the spin susceptibility would keep the full value 
of the normal-state Pauli spin susceptibility for magnetic 
fields along the c-axis. Therefore, no change of phase 
would be induced by a magnetic field. For the interme- 
diate situation, if spin singlet and spin triplet pairing are 
of comparable strength an intriguing mixing of the two 
channels is possible leading to a complex phase structure 
involving spontaneous time reversal symmetry breaking. 
This situation will be discussed elsewhere. Note, how- 
ever, that the pairing interactions in the spin triplet chan- 
nel are expected to be weak in the bulk/superlattice of 
CeCoIns and high-Tc cuprates. Even then, weak interac- 
tions in the spin triplet channel gives rise to the uniform 
spin triplet Cooper pairs in the PDW state. 

Second, we focused on the Pauli-limited superconduc- 
tor with large Maki parameter, neglecting orbital depair- 
ing effects, in particular, the aspect of a mixed phase in- 
cluding flux lines. Although orbital depairing may play 
quantitatively important roles for the field along c-axis, 
the PDW state survives moderate orbital depairing ef- 



fects (moderate Maki parameters) because the spatial 
modulation of vortex state along the a6-plane is not af- 
fecting the TT-phase shift along c-axis profoundly. The 
situation is quite different for the helical superconducting 
state in Rashba-type non-centrosymmetric superconduc- 
tors, which is obscured even by the weak orbital depairing 
effect.?? 

Third, the PDW state is stabilized not only in the bl- 
and tri-layers but also in systems with more layers. In 
such systems the layers towards the center are subject 
to paramagnetic limiting and it is expected that within 
the BCS state we encounter a progressive demolishing of 
superconductivity in these layers upon increasing field. 
Eventually, at high enough magnetic field the supercon- 
ducting system behaves like the bi-layer system in the 
PDW state. Due to the fact that these outer layers are 
well separated the condition for the appearance of the 
PDW state is less stringent than in the real bi-layer sys- 
tem. These views have been confirmed by us numerically 
for the 4-layer system. 

In conclusion, we have studied the superconducting 
state in multilayer systems including layer dependent an- 
tisymmetric spin-orbit coupling due to the coordination 
of the layers. We found that in a magnetic field along 
the c-axis spin-orbit coupling influences the phase of the 
order parameters of different layers. We find a first or- 
der phase transition separating two states different by 
symmetry: the high-field BCS state which is even un- 
der reflection at the center of the multilayer system and 
the PDW state which is odd under reflection. The latter 
phase allows for a considerable enhancement of the upper 
critical field and generates the uniform order parameter 
for the spin triplet superconductivity. This phenomenol- 
ogy relies on the fact that orbital depairing is weak, i.e. 
the coherence length of the superconductor is very short, 
and that the dominant pairing is in a spin singlet channel. 
This may apply to the artificially grown multilayer ver- 
sion of CeCoIus and some multilayered high-Tc cuprates. 
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Appendix: CASE OF THE D-WAVE 
SUPERCONDUCTIVITY 

Here we show the phase diagram for d-wave supercon- 
ductivity using a linearized gap equation. We define the 
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superconducting susceptibility Xmim2il) following, 



Jo 



Sq„,(r)Sq„,(0)), 



(A.l) 



where q = (q,iri„), r2„ = 27rn//3 with an integer 
n is a boson Matsubara frequency, fi = l/k-oT, and 



is the creation operator of spin-singlet Cooper pairs 
with a center-of-mass momentum q on the layer m, and 
(j){k) = cos fcx — cos fcj, shows the momentum dependence 

I 



of the order parameter with d-wave symmetry. Tak- 
ing into account the attractive interaction in the d-wave 
channel Vd, we obtain the following expression. 



1 - VdXtiq) 



(A.2) 



where x^'^iq) and yt^iq) are the M x M matrix repre- 
sentation for the reducible and irreducible superconduct- 
ing susceptibility, respectively {M is the number of lay- 
ers). The bare superconducting susceptibility Xo^{q) is 
obtained as 



Xor, 



k I 



~Gl\,,„,^{k + q,iui)G\l^„^^{-k,inn - iuji)<j){k)(j)*{-k - g)], 



(A.3) 



where uji — Tr{2l + is a fermion Matsubara frequency 
and Gf^^j^^ {k, iui) is the non- interacting Green function. 

The superconducting critical temperature is deter- 
mined by the divergence of the superconducting suscep- 
tibility x"'=(0,0), namely i - VdXa{ld,Q) = (Thouless 
criterion). 

In Fig. [SJ we show the phase diagram for bi- and tri- 
layer systems. Although we cannot determine the first 



order phase transition line between the BCS and PDW 
states within the linearized equation, the superconduct- 
ing transition lines in Fig. IHJa) and Fig. IHJb) resemble 
those in Fig.[2ljd) and Fig.[3j respectively. Hence, our re- 
sults in this paper are qualitatively valid for the d-wave 
superconducting state. 



1 P. Fulde and R. A. Ferrell, Pliys. Rev. 135, A550 (1964). 

^ A. I. Larkin and Yu. N. Ovchinnikov, Zh. Eksp. Teor. Fiz. 
47, 1136 (1964) [Sov. Phys. JETP 20, 762 (1965)]. 

^ H. A. Radovan, N. A. Fortune, T. P. Murphy, S. T. Han- 
nahs, E. C. Palm, S. W. Tozer, and D. Hall, Nature (Lon- 
don) 425, 51 (2003). 

* A. Bianchi, R. Movshovich, C. Capan, P. G. Pagliuso, and 
J. L. Sarrao, Phys. Rev. Lett. 91, 187004 (2003). 

^ Y.-a. Liao, A. S. C. Rittner, T. Paprotta, W. Li, G. B. 

Partridge, R. G. Hulet, S. K. Baur, and E. J. Mueller, 

Nature (London) 467, 567 (2010). 
^ R. Casalbuoni and G. Nardulli, Rev. Mod. Phys. 76, 263 

(2004). 

"Non-centrosymmetric Superconductors", edited by E. 
Bauer and M. Sigrist (Springer- Verlag, 2012). 

* Y. Matsuda and H. Shimahara, J. Phys. Soc. Jpn. 76, 
051005 (2007). 

^ K. Kumagai, H. Shishido, T. Shibauchi, and Y. Matsuda, 

Phys. Rev. Lett. 106, 137004 (2011). 
^° Y. Mizukami, H. Shishido, T. Shibauchi, M. Shimozawa, 
S. Yasumoto, D. Watanabe, M. Yamashita, H. Ikeda, 
T. Terashima, H. Kontani, and Y. Matsuda, Nature Phys. 
7, 849 (2011). 

D. Maruyama, M. Sigrist, and Y. Yanase, J. Phys. Soc. 
Jpn. 81, 034702 (2012). 
^2 Y. Yanase, J. Phys. Soc. Jpn. 79, 084701 (2010). 

S. K. Goh, Y. Mizukami, H. Shishido, D. Watanabe, S. 
Yasumoto, M. Shimozawa, M. Yamashita, T. Terashima, 



Y. Yanase, T. Shibauchi, A. L Buzdin, and Y. Matsuda, 
arXiv: 1207.4889. 
" J-H. She and A. V. Balatsky, Phys. Rev. Lett. 109, 077002 
(2012). 

Y. Yanase and M. Sigrist, J. Phys. Soc. Jpn. 77, 124711 
(2008). 

^® D. F. Agterberg, M. Sigrist, and H. Tsunetsugu, Phys. 
Rev. Lett. 102, 207004 (2009). 

A. Aperis, G. Varelogiannis, and P. B. Littlewood, Phys. 
Rev. Lett. 104, 216403 (2010). 

For a review, see A. L Buzdin, Rev. Mod. Phys. 77, 935 
(2005). 

H. Mukuda, M. Abe, Y. Araki, Y. Kitaoka, K. Tokiwa, 
T. Watanabe, A. lyo, H. Kito, and Y. Tanaka, Phys. Rev. 
Lett. 96, 087001 (2006). 
^" S. Shimizu, S. Tabata, H. Mukuda, Y. Kitaoka, P. M. Shi- 
rage, H. Kito, and A. lyo, J. Phys. Soc. Jpn. 80, 043706 
(2011). 

21 J. L. O'Brien, H. Nakagawa, A. S. Dzurak, R. G. Clark, 

B. E. Kane, N. E. Lumpkin, N. Miura, E. E. Mitchell, 
J. D. Goettee, J. S. Brooks, D. G. Rickel, and R. P. Star- 
rett, Phys. Rev. B 61, 1584 (2000). 

22 Y. Matsunaga, N. Hiasa, and R. Ikeda, Phys. Rev. B 78, 
220508 (2008). 



7 



Normal 




(b) 



X 

CD 





iNormai 


Superconductivity 





0.2 



0.4 0.6 

T/T, 



1.2 



CO 



FIG. 6. T-H phase diagram for (a) bi- and (b) tri-layer system, with attractive interaction in the d-wave symmetry, V{k, k') — 
— Vd{coskx — cos fcy)(cos fe^ — cosk'y). We assume Vd/t = 2.2 and a/t± = 3. The critical temperature is kBTco/t = 0.0201 for 
bi-layer system and kBTco/t = 0.0202 for tri-layer system for a = 0. Other parameters are the same as in Fig. (2] The points 
show the end point of first order transition line between the BCS state to the PDW state. 



